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We investigate the influence of AC driving fields on the coherence properties of one- and two-qubit gate 
operations. In both cases, we find that for suitable driving parameters, the gate purity improves significantly. 
A mapping of the time-dependent system-bath model to an effective static model provides analytical results. 
The resulting purity loss compares favorably with numerical results. 



1 Introduction 

The experimental realization of one-qubit gates in solid state setups [1-3] and two-qubit gates in ion 
traps [4, 5] and Josephson junctions [6] has demonstrated that these systems provide remarkable coherence 
properties although the goal of 10~ 5 errors per gate operation [7] has not yet been accomplished experi- 
mentally. The unavoidable coupling to external degrees of freedom and the thereby caused decoherence 
still presents a main obstacle for the realization of a quantum computer. Several proposals to overcome 
the ensuing decoherence have been put forward, such as the use of decoherence free subspaces [8-12], 
coherence-preserving qubits [13], quantum Zeno subspaces [14], optimized pulse sequences [15, 16], dy- 
namical decoupling [17-21], and coherent destruction of tunneling [22,23]. Theoretical studies of deco- 
herence of two-level systems have been extended to gate operations in the presence of an environment in 
Refs. [24-29]. 

A variety of suggestions towards coherence stabilization relies on the influence of external fields. One 
of the most prominent examples is the application of a sequence of 7r-pulses that flip the sign of the qubit- 
bath coupling operator resulting in a so-called dynamical decoupling (DD) of the qubit from the bath 
[17-21,23]. A drawback of this scheme is the fact that it eliminates only noise sources with a frequency 
below the repetition rate of the pulses. This clearly causes practical limitations. However, these limitations 
may be circumvented by using a related scheme based on continuous-wave driving, i.e. one with a harmonic 
time-dependence, which allows higher driving frequencies. 

A different proposal for coherence stabilization is to employ the physics of the so-called coherent de- 
struction of tunneling (CDT). CDT has originally been discovered in the context of tunneling in a driven 
bistable potential [30-34]. There, it has been found that a particle which is initially in the, say, left well of 
a symmetric bistable potential, can be prevented from tunneling by the purely coherent influence of an os- 
cillating driving field. This effect is stable against dissipation in the sense that the AC field also decelerates 
the dissipative transitions from the left to the right well [35-37]. 

A less frequently studied problem is the extension of coherence stabilization protocols to systems that 
consist of two or more interacting qubits. Such an interaction is essential for two-qubit operations which 
represent an indispensable part of all quantum algorithms [38^-0]. For a CNOT gate based on an isotropic 
Heisenberg interaction [24,41], it has been proposed to stabilize coherence by applying a control field to 
one of the qubits and thereby obtain an effective Ising interaction which is less sensitive to the influence of 
environmental degrees of freedom [42]. This scheme possesses the beneficial properties that (i) it involves 
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only intermediately large driving frequencies that can lie well below the bath cutoff and (ii) it does not 
increase the gate operation time. Moreover, since the driving field couples to the same coordinate as the 
quantum noise, this coherence stabilization is distinctly different from the recently measured dynamical 
decoupling of a spin pair from surrounding spin pairs [43]. 

In this work, we extend our previous studies on coherence stabilization of one- [23] and two-qubit 
operations [42] and, moreover, detail some technical aspects. In Section|2] we introduce a model for two 
qubits coupled to a heat bath and derive a Bloch-Redfield master equation to describe quantum dissipation 
and decoherence. This formalism is applied in Sections[3]and|4] respectively, to the single qubit dynamics 
and to a two-qubit gate operation. In each case, we derive within a rotating-wave approximation (RWA) 
analytical results for the decay of the gate purity. Thereby, the proper treatment of the qubit-bath coupling 
is of crucial importance. The computation of averages for the ensemble of all pure states is deferred to the 
appendix. 



2 Quantum gate with bit-flip noise 

We consider a pair of qubits described by the Hamiltonian [24, 38^-1] 

ffqubits = \^2 ( A J' CT J + e ^ X ) + J(? i ' ^2, (1) 

3=1,2 

with a qubit-qubit coupling of the Heisenberg type, where j labels the qubits. In order to construct a 
quantum gate, the tunnel splittings Aj, the biases Ej, and the qubit-qubit coupling J have to be controllable 
in the sense that they can be turned off and that their signs can be changed. 
The bit-flip noise is specified by the system-bath Hamiltonian [44,45] 

H = -ffqubits + -ffcoupl + -ffbath (2) 

where 

-ffcoupl = - a 3 hc »( a ]u + a 3v) ( 3 ) 

3=1,2 u 

denotes the coupling of qubit j to a bath of harmonic oscillators with frequencies uj v described by the 
Hamiltonian Hbath = J2ju f^i/dj^o-ju and the spectral density I(u>) = it cl5(u> — cj v ). This coupling 
of each qubit to an individual bath represents a proper model for sufficiently distant qubits [29]. Within the 
present work, we consider the so-called ohmic spectral density 

I(lu) = (4) 

with the dimensionless coupling strength a and the cutoff frequency lo c . In order to complete the model, we 
specify the initial condition of the Feynman- Vernon type, i.e., initially, the bath is in thermal equilibrium 
and uncorrelated with the system, ptot(io) = p(to) ® -Rbath,cq, where p denotes the reduced density 
operator of the two qubits and i?bath,cq oc exp(— -ffbath/fcB? 1 ) is the canonical ensemble of the bath. 



2. 1 Bloch-Redfield master equation 

If ak-g,T is smaller than the typical system energy and if the dissipation strength is sufficiently small, 
a 1, the dissipative system dynamics is well described within a Bloch-Redfield approach, which is 
also referred to as Born-Markov approach. There, one starts from the Liouville-von Neumann equation 
iSptot = [H, ptot] for the total density operator and obtains by standard techniques the master equation [46] 

P = - ^[H quhits , p] - [Q {t), p]] - £>J, {Pj(t),p}} (5) 

3 3 

= - ^[#qubit SJ />]-A(t)p (6) 
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with the anti-commutator {A, B} = AB + BA and 

-i poo fOO 

Q j (t) = — dr / dwS(w)cos(wT)or?(t-T,t), (7) 

■t poo poo 

P j (t) = — dr du)I(u))sin(uT)a?(t-T,t). (8) 

47r Jo Jo 

Thus, the influence of the bath is determined by the Heisenberg operators of the system, the spectral density 
I(oj) of the heat baths, and the Fourier transformed 

S(lo) = I(lu) coth(hLo/2k B T) (9) 

of the symmetrically-ordered equilibrium auto-correlation function ^({^ (r), £j (0)}) eq of the collective 
bath coordinate £j — c "i a ju + a jv)- The notation X(t, t') is a shorthand for the Heisenberg operator 
U'(t, t')XU(t, t') with U being the propagator of the coherent system dynamics. Note that S(ui) and I(uj) 
are independent of j due to the assumption of two identical environments. We emphasize that the particular 
form of the master equation is valid also for an explicitly time-dependent qubit Hamiltonian. 



2.2 Purity decay 

The heat baths, whose influence is described by the second and third term of the master equation 0, lead to 
decoherence, i.e., the evolution from a pure state to an incoherent mixture. This process can be measured 
by the decay of the purity tr(p 2 ) from the ideal value 1. The gate purity (frequently also referred to as 
"purity") V(t) = tr(p 2 (t)), which characterizes the gate independently of the specific input, results from 
the ensemble average over all pure initial states [47]. For weak dissipation, the purity is determined by its 
decay rate at initial time, 

r =-^lnP(t)|^ = 2te^= ] ^^Str(aJ0 i ). (10) 

In order to obtain the last expression, we have used the cyclic property of the trace and performed the 
ensemble average over all pure states as described in AppendixIXl Here, N is the dimension of the system 
Hilbert space. The purity loss rate V represents a figure of merit for the coherence of the quantum gate — 
ideally, it vanishes. Interestingly enough, only the second term on the right-hand side of the master equation 
(|3J contributes to the purity decay. This relates to the interpretation that the first term of the master equation 
is responsible for the coherent dynamics, while the second and third term correspond to decoherence and 
relaxation, respectively. 

The average over all pure states, which underlies the rate ( llOl l, may differ from the average over a 
discrete set of initial states by a factor of the order unity, due to the particular choice of a discrete ensemble: 
There, one commonly chooses ensembles which are unsymmetric on the Bloch sphere [27, 28] or which 
do not include entangled states [27-29]. 



2.3 Numerical solution 

The purity decay rate dlOi by construction accounts only for the behavior at initial time t = 0. Thus, 
for a more complete picture, it is desirable to have the exact numerical solution of the master equation 
(JSJl at hand. For studying the influence of an external AC field, such a solution must properly capture 
the case of a T-periodic system Hamiltonian. An efficient scheme for that purpose is a modified Bloch- 
Redfield formalism whose cornerstone is a decomposition into the Floquet basis of the driven system 
[46]: According to the Floquet theorem, the Schrodinger equation of a driven quantum system with a 
Hamiltonian of the form H(t) = H(t + T) possesses a complete set of solutions of the form \ip a (t)) = 
exp(—ie a t/h)\<f> a (t)). The so-called Floquet states \<j> a (i)) obey the time-periodicity of the Hamiltonian 
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Fig. 1 Bloch sphere representation of the single-qubit Hamiltonian 1151 and the bath coupling studied in 
Sec.[3] 



and e a denotes the so-called quasienergy. The Floquet states are elements of an Hilbert space extended by 
a T-periodic time coordinate and are computed from the eigenvalue equation 

(H(t)-ih^)\<t> a (t))=e a \<f> a (t))- (ID 
In the Floquet basis { 1 4> a (t)) }, the master equation Q assumes the form 

Pa(3 = — jX e a — £p)PaP ~ ^ & a p,a'P'(t) Pa' 0' ■ (12) 

a'P' 

Besides computational advantages, using the Floquet basis has the benefit that it implies the numerically 
exact treatment of the coherent system dynamics. Thereby, one avoids artefacts like the violation of equi- 
librium properties in the undriven limit [48-52]. 

For weak dissipation, we can replace within a rotating-wave approximation A(<) by its time aver- 
age [46]. Finally, we integrate the master equation to obtain the dissipative propagator W a p,a'P' which 
provides the final state p ou t, a p = Y^a'P' W a p,a'P'pin,a'P' an d, thus, allows one to evaluate all quanti- 
ties of interest. In particular, it is possible to compute for a pure initial state \tp) = J2 a c a\4>a)i i- e - f° r 
Pin, ap — c a c*p, the final state p out which possesses the purity 

tr Pout = Pout, aP Pout, Pa = ^ ~WaP,a'P' ~Wpa,P"a" C a > C*p,Cp» C*« . (13) 

ap aPa' P' a" P" 

In order to average over all pure initial states, we employ Eq. d48l derived in the Appendix to obtain the 
gate purity 

Pout = N(N s (w a p, a <p>Wp a ,p, a , +W a p, a ' a -Wp a ,p>p,y (14) 

( ' aPa'P' 

We emphasize that this result is independent of the particular choice of the basis. 

3 Coherence stabilization for single qubits 

During the stage of single qubit operations, i.e., for J = 0, both qubits together with the respective bath 
evolve independently of each other. Thus, it is sufficient to focus on qubit 1 with the Hilbert space dimen- 
sion N = 2. We restrict ourselves to an operation with e% = and the fixed tunnel splitting Ai > 0, i.e. to 
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the Hamiltonian 

Ai 

-Hqubitl = _ 2 _<J f (1^) 

This situation is characterized by the fact that the bath couples to a system operator which is different from 
the qubit Hamiltonian. In the Bloch sphere representation sketched in Fig. [2 the qubit Hamiltonian dl5> 
and the bath coupling correspond to orthogonal vectors. We start out by bringing the master equation (0 
into a more explicit form by inserting the Heisenberg operator 

of(t- r,t) = of cos ^j— +of sin— jj^, (16) 

which is readily derived from its definition together with the qubit Hamiltonian Q. Performing the integra- 
tion over t and ui and neglecting renormalization effects, which are small provided that a ln(oj c / Ai) <C 1, 
yields for Ai <C uo c , the Markovian master equation 

P = ~ ~ [ffqubiti, P] ~ [of , [of , p]} + i^i [of , {o\, p}} . (17) 

Then the purity decay dlOt is readily evaluated to read 

r = ls(A/H), (18) 
6 

where the subscript "0" refers to the absence of any AC field. 

A possible coupling to an external driving field might have any "direction" n on the Bloch sphere (cf. 
Fig.[0, i.e. be proportional to n ■ a%. Herein, we restrict ourselves to the cases parallel to the bath coupling 
and parallel to the static Hamiltonian, i.e., to a field that couples to of and of, respectively. 

3.1 Dynamical decoupling by harmonic driving 
The first option is to act on the qubit with a driving of the form 

A 

#11 = -o\ cos(Ot) (19) 

"parallel" to the static Hamiltonian. This relates to a recently proposed mechanism for coherence stabi- 
lization, namely the so-called dynamical decoupling (DD) [17-21]. This scheme employs sequences of 
7r-pulses that flip the sign of the operator of which couples the qubit to the bath. The basic idea dates back 
to the suppression of spin diffusion in nuclear magnetic resonance experiments [53, 54] and by now is an 
established technique in that area [55]. In the present case where the bath couples to the operator of [cf. 
Eq. (|3J], such a transformation is e.g. induced by the Hamiltonian Hljr .of for a pulse duration tv/u)r. Since 
the corresponding propagator is a function of the qubit Hamiltonian, the coherent dynamics is not altered. 
Besides the prospective benefits of such a control scheme, there is also a number of possible drawbacks that 
the application of it -pulses might cause: For a driven system, there is always the possibility of unwanted 
off -resonant transitions [56], especially in the case of ideal rectangular pulses. A more practical limitation 
is the fact that only noise with frequencies below the pulse repetition rate can be eliminated in such a way. 
These disadvantages can be overcome partially by applying a continuous wave version of the dynamical 
decoupling scheme, i.e. a driving of the form dl9l for which the available frequency range is larger. 

For the computation of the coherence properties, we use the fact that Hu (t) commutes with the static 
qubit Hamiltonian d!5l > and, consequently, the propagator for the driven qubit can be computed exactly 
reading 

U(t,t') = exp (-' 1 2^[ s H nt ) - sin(QOK ) ex P (-^Ai erf (i - t'U . (20) 
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Fig. 2 Decoherence reduction by dynamical decoupling, 77dd 5 as a function of the driving frequency for various 
temperatures. The cutoff frequency is u c = 500Ai//i, A/hQ = 2.4, and the dissipation strength is a — 0.01. The 
horizontal line marks the value 1 below decoherence is lower than in the static case. 



We have written the propagator in a form that is suitable for simplifying the master equation (|5}- Inserting 
this into the expression Q results for Ai <C Jl in the effective coupling operator 

Qdd = \ ^(A/hQ)S(A 1 /h) + 2 Jl (A/hQ) S(nJl)^ of. (21) 

In order to derive this expression, we have decomposed the exponentials of the trigonometric functions 
into a Fourier series using the identity exp[ixsin(f2t)] = J2k Jk(x) exp(ifcSlt), where Jk is the fcth order 
Bessel function of the first kind [57]. The effective coupling operator Qdd is proportional to a'f and, thus, 
the master equation is again of the form Q. The only difference is that the dissipative terms have acquired 
the prefactor 

which modifies the purity decay rate accordingly, Tdd = ^ddTo- Equation \22\ allows for the interpreta- 
tion that now the decoherence rate depends on the spectral density of the bath at multiples of the driving 
frequency which may be larger than the cutoff frequency ui c . The 7r-pulses applied in the original ver- 
sion [17] of dynamical decoupling, correspond for a continuous driving to a field amplitude that is adjusted 
such that A/Ml equals the first zero of the Bessel function Jo, i.e. it assumes a value 2.404825 . . .. Then 
only the sum in Eq. J22i contributes to the decoherence rate Tdd ■ If now the driving frequency is larger 
than the cutoff of the spectral density, SI > uj c , decoherence is considerably reduced: For low temperatures, 
k B T -c Ai, the hyperbolic tangent in the factor 1221 becomes unity and each contribution is weighted by 
a possibly large factor nhfl/Ai. In the high-temperature limit k B T ^> Ml, we use tanh(x) x and find 
that the dependence of the prefactor on nfl cancels. This means that the dynamical decoupling scheme 
is especially useful for high temperatures. The physical reason for this is that the driving shifts the qubit 
dynamics towards high frequencies where the thermal occupation of the bath modes is negligible. 

Figure|2]compares the coherence stabilization t^dd as a function of the driving frequency for A/Ml = 
2.4, i.e. close to a zero of the Bessel function Jo. It reveals that for driving frequencies well below the 
cutoff, the driving rather spoils the coherence. This improves with increasing driving frequency and, 
finally, for a high-frequency driving, ?7dd becomes much smaller than unity corresponding to a significant 
coherence stabilization. The data demonstrate the particular usefulness of dynamical decoupling at high 
temperatures. 
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3.2 Coherent destruction of tunneling 
Our second example under consideration is a driving field 

H ± (t)=f(t)af, (23) 

where the energy / (t) is a 27r/ -periodic function of time with zero mean. The field couples to the qubit by 
the same operator erf as the bath, i.e., "perpendicular" to the static Hamiltonian dl5> . Thus, it commutes 
with the qubit-bath coupling but not with the static Hamiltonian. Such a time-dependent field causes 
already interesting effects for the coherent qubit dynamics that we will briefly review before discussing 
decoherence. 

For that purpose, we derive within a rotating-wave approximation (RWA) analytical expressions for 
both the coherent propagator U(t, t') and the purity decay JlOi . We start out by transforming the total 
Hamiltonian into a rotating frame with respect to the driving via the unitary transformation 

CacW = e-W*)"*, <f>(t) = ^j\t'f(tf). (24) 
This yields the likewise 27r/Sl -periodic interaction-picture Hamiltonian 

ff q ubiti(*) = t/lcWtfqubitif/AcW (25) 
Ai f , ( -1 \ , v ... (A 



= T i 1 cos [m sin(m) ) + ^ sin [m sin(m) ) ) (26) 

and the S-matrix S(t, t') = U AC {t)U{t,t')UAc{t')- The corresponding Schrodinger equation cannot be 
integrated exactly since if qu biti(i) does not commute with itself at different times and, thus, time-ordering 
has to be taken into account. We restrict ourselves to an approximate solution and neglect corrections 
of the order A^. Within this approximation, the propagator is simply given by the exponential of the 
integral of the time-dependent interaction-picture Hamiltonian. This is equivalent to replacing \26\ by its 
time-average 

-ffqubitl = (-Hqubitl (*))27r/n = — (27) 

where (. . .^/n denotes the time-average over the driving period. This RWA approximation to the driven 
qubit Hamiltonian is of the same form as the original static Hamiltonian Q, but with the tunneling matrix 
element being renormalized according to 

Ai A off = (cos[20(i)]) 2x/o A 1 , (28) 

Consequently, we find S(t, t') = exp{— iH qu biti(t ~ t')/h}, such that within RWA, the entire propagator 
for the qubit in the Schrodinger picture reads 

U(t,t') — e~ i< ^ t ) CT i e ~iH^hiti(t~t')/h & i(j>(t')af ^ ^29) 



Of particular interest are now driving parameters for which the effective tunnel splitting (12 8 1 and, thus, 
^qubiti vanish. Then, the one -period propagator U(t+T, t) becomes the identity [recall that Uac is 27r/f2- 
periodic and, thus, [/ac(27t/^) = Uac(0) = 1]. This implies that the long-time dynamics is suppressed. 
The dynamics within the driving period requires a closer look at the 27r/Sl-periodic contribution UAc{t)'- 
For an initial preparation in an eigenstate of af, it provides only a global phase, such that the dynamics 
as a whole is suppressed also within the driving period. This effect of suppressing the time-evolution 
by the purely coherent influence of an external field has been investigated first in the context of driven 
tunneling [30, 31] and is named "coherent destruction of tunneling" (CDT). Therefore, we will refer to a 



8 



S. Kohler and P. Hanggi: Improving the purity of one- and two-qubit gates by AC fields 



driving of the form d23i as "CDT driving" despite the fact that we also consider working points at which 
the coherent dynamics is not completely suppressed. Note that for a preparation other than an eigenstate 
of erf, the periodic propagator U^cif] will still cause a non-trivial dynamics within the driving period. 

Let us now turn to the influence of the CDT driving (I23l > on quantum dissipation and decoherence. For 
that purpose, we have to evaluate the operator Qj contained in the master equation (|5jl. Inserting the RWA 
propagator J29t into 0, we obtain after some algebra the result 

Qcdt = ^S(A cS /h)af. (30) 

Thus, the master equation again is of the same form as in the undriven case, Eq. l|5}, but the generator of 
the dissipative dynamics is modified by the factor 

S(A cS /h) coth(A cff /2fc B T) , „., ffl nu 
?7cdt = = coth(Al/2fcBr) (cos[2#)]) 27r/n . (31) 

Consequently, the purity decay becomes Tcdt = ^cdtTo- Since the spectral density S(ui) increases 
monotonically with the frequency u> and, moreover, the Bessel function Jq(x) < 1, the CDT driving — in 
clear contrast to the dynamical decoupling — never enhances dissipation and decoherence. 

In the high-temperature limit k-&T ^> A e fi, we employ the approximation coth(x) » 1/x which 
implies that for an ohmic bath, S(A/h) w AnakBT/H is independent of the tunnel splitting. Consequently, 
Vcdt ~ 1, i-e., the purity decay is essentially unchanged. 

In the opposite limit of low temperatures, k^T <C A c ff, the argument of the hyperbolic cotangent 
is large such that coth(x) w 1. Then ?7cdt = A c ff/Ai = (cos[20(f)])27rn < 1- This reduction of 
decoherence is brought about by the fact that the driving d23l > decelerates the long time dynamics of the 
qubit. Thereby, the frequencies which are relevant for the decoherence are shifted to a range where the 
spectral density of the bath is lower. Consequently, the influence of the bath is diminished. 

An important special case is that of a harmonically time-dependent driving field, f(t) = ^Acos(ilt). 
Then, the time-average in the effective matrix element J28t can be expressed by a zeroth-order Bessel 
function of the first kind such that A e g = Jq(A/MY)Ai. Consequently, we find that at sufficiently low 
temperatures decoherence is reduced by a factor ?7cdt = Ja(A/Ml). These results imply that for an ohmic 
bath, the coherent dynamics is slowed down by the same factor as the decoherence, cf. Eqs. (12 8 1 and J3 1 1 in 
the low-temperature limit. Thus, if for a specific application, the figure of merit is the number of coherent 
oscillations, the present coherence stabilization scheme may therefore not prove very useful. 



4 Coherence stabilization for a CNOT gate 

Logic operations in both quantum and classical computers require that the time evolution of a qubit (or a 
bit, respectively) depends on the state of another qubit. The physics behind such processes is a non-linear 
interaction between two qubits which in many implementations is of the Heisenberg type [24,41] 

-^Heisenberg = J <?1 ' <?2 (32) 

as assumed in our model Hamiltonian Q. In the context of decoherence, the question arises whether 
it is possible to stabilize coherence by a proper AC field while maintaining the desired non-linear time 
evolution. 

The qubit Hamiltonian Q together with the bath coupling (0, allows one to implement the CNOT 
operation sketched in Figure [3J> [24, 27, 38^-0, 58]. Apart from single qubit operations, it consists of the 
propagator 

Un(f) = exp(-i<p<?i • ct 2 ) (33) 

for the Heisenberg qubit-qubit interaction (1321 which in total acts for a time t j = ttH/AJ such that ip = 
7r/4. For single qubit operations, we have discussed in Sec.f3]that pulse sequences [17-21] and harmonic 
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Fig. 3 Realizations of a CNOT operation for different types of qubit-qubit interaction [38-40]. (a) Heisen- 
berg interaction providing the transformation 1331 and (b) Ising coupling along the x-direction correspond- 
ing to I42i . The symbol [H] denotes the Hadamard operation and [4>] n the rotation of the respective qubit 
around the axis n by an angle cj>. 



driving fields [23] can suppress decoherence. Therefore, we focus here on decoherence during the stage 
of the qubit-qubit interaction and, thus, take as a working hypothesis that the coherence of one-qubit 
operations can be stabilized ideally. Then the remaining decoherence takes place during the qubit-qubit 
interaction time t j. 



4. 1 Heisenberg vs. Ising coupling 

Before discussing the influence of a driving field, we provide for later reference the results for the purity 
decay in the static situation. Therefore, we need to evaluate tr (a'^Qj) where 

-| /*00 /"OO 

Q ] = — dt dujS(cu)cos(LUT)e iHK ^ b ^ T/h a^e~ iHK ^ b ^ T/h . (34) 
4tt Jo Jo 

This calculation is most conveniently performed in the basis of the total (pseudo) spin L = + a 2) 
because i?Hciscnbcrg = J(2L 2 — 3), where L 2 possesses the eigenvalues £(£ + 1), £ = 0, 1. For £ = 0, the 
corresponding eigenstate is the singlet state (|01) — 1 10} ) / a/2 with energy eo = —3 J, while for I = 1, one 
finds the triplet |00), (|01) + |10))/\/2, 1 1 1 > with ei = J. After evaluating the matrix elements of crj and 
Qj, we finally arrive at the purity decay rate 

THciscnbcrg = ^{5(0) + 5(4 J/fi)} (35) 

which is sketched in Fig. [4] In the derivation, we have ignored Lamb shifts and defined 

47T 

5(0) = lim S(uj) = —ak B T. (36) 

In particular, we find that for low temperatures, U-qT < J, decoherence is dominated by the term 5(4 J/ H) 
such that r ps 167raJ/5fi. This part reflects the influence of the so-called quantum noise which is 
temperature-independent and, thus, cannot be reduced by cooling the environment. 

It is intriguing to compare the result M5\ with the purity decay of a similar system, namely a pair of 
qubits interacting with the Ising interaction 

Rising = Jo\a x 2 . (37) 
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Fig. 4 Purity loss (1 — V) of two qubits interacting for a time t j = nh/4J comparing Heisenberg (solid) 
and Ising (dashed) coupling. The dimensionless dissipation strength is a — 0.01/27T. 



Such an interaction can also be employed for the construction of quantum gates like the CNOT gate in 
Fig.[3p. (Note, however, that for several physical realizations of a quantum computer [24,41], the natural 
interaction is of the Heisenberg type). The main difference to the Heisenberg interaction ( I32l i. is the fact 
that i/ising commutes with the bath coupling Q, i.e., the bath couples to a good quantum number. Carrying 
out the same calculation as above yields the purity decay 

ri S i„ g = ^5(0) = liraknT (38) 
5 5 

which exhibits a significantly different low-temperature behavior: Instead of saturating, it remains propor- 
tional to the temperature, cf. Fig.|4] 



4.2 Coherence stabilization by an AC field 

In Sections 13 ,2l and l4~T1 we presented two results that lead to the central idea for coherence stabilization of 
a qubit pair under the influence of a Heisenberg coupling: First, an AC field can suppress the dynamics 
"transverse" to the driving, i.e., it can effectively eliminate the parts of the Hamiltonian that depend on spin 
matrices other than erf . Second, qubits with Ising interaction are less sensitive to decoherence than qubits 
with Heisenberg interaction. Thus, the question arises whether one can act with an AC field on the system 
in such a way that precisely the part of the Hamiltonian d32t that causes the quantum noise becomes 
suppressed. This is indeed the case and can be performed by driving qubit 1 with the AC field i23\ . Note 
that qubit 2 remains undriven. 

For the computation of the coherence properties, we proceed as in Section l3~2l We first derive an 
effective static qubit Hamiltonian by transforming the time-dependent Hamiltonian i?Heisenberg + H± (t) 
via (1241 into rotating frame and subsequently replace it by its time average. After some algebra along the 
lines of Section l3~2l we obtain the Hamiltonian 

-ffqubits = (J~ Jjjofof +J±Vl-&2l ( 39 ) 

where the constant 



Jj_ = J(cos[2^(t)]) 27r/Q (40) 

denotes an effective interaction "transverse" to the driving and (. . .^wn the time average over the driving 
period. Consequently, we find the S-matrix S(t, t') = exp{— i-ff qu bits(i — t')/K\, such that the propagator 
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of the driven system again assumes the form d29i with i/ qu biti replaced by -ff qu bits- Having this propagator 
at hand, we are in the position to derive explicit expressions for the operators erj(i — r, t) and Qj, Again, 

the calculation is conveniently done in the basis of the total spin L and L x which, owing to the relation 
erf erf — i (erf + erf ) 2 — 1, is an eigenbasis of the Hamiltonian 09l >. We evaluate the purity decay rate dlOt 
in this basis and finally obtain 

Hciscnbcrg, driven — 

-{S(0)+S(4J ± /h)}, (41) 

i.e., the result d35l but with J replaced by J±. For f(t) = 0, we find J± = J such that the static result is 
reproduced; otherwise, the inequality \J±\ < J holds and, thus, the bath correlation function S in Eq. J41i 
has to be evaluated at a lower frequency. For an ohmic or a super-ohmic bath, S(ui) is a monotonously 
increasing function and, consequently, the AC field reduces purity decay (unless J > cj cu t ff)- 

The purity decay assumes its minimum for J± = 0. This condition marks the working points on which 
we shall focus henceforth. For an ohmic spectral density I(uj) = I-koluj, the purity decay at the working 
points becomes T = |5(0) = Snak^T /b. This value has to be compared to the purity decay in the 
absence of driving: An analysis reveals that for k-gT > J, decoherence is essentially driving independent. 
By contrast for low temperatures, k-gT < J, the driving reduces the decoherence rate by a factor k-gT /2 J. 
This low-temperature behavior results from the fact that for J± = 0, the effective Hamiltonian d39l is iden- 
tical with the Ising Hamiltonian J37i and, thus, commutes with the qubit-bath coupling operators crj. Note 
that the latter are not affected by the transformation i24l . This means that the driving modifies the effective 
qubit Hamiltonian such that the bath acts as pure phase noise whose influence is proportional to the tem- 
perature. The fact that the driving field couples to the same coordinate as the bath distinguishes the present 
coherence stabilization from dynamical decoupling. In that respect, the present scheme is complementary 
to coherence-preserving qubits [13], for which heating errors are the only source of decoherence. 

For a rectangular driving for which f(t) switches between the values ±A/2, the condition J± = 
is equivalent to A = Ml and corresponds to two 7r-pulses per period. For a harmonic driving, f(t) = 
Acos(ilt)/2, one obtains J± = J Jq(A/MV), where Jo denotes the zeroth-order Bessel function of the 
first kind. Then, at the working points J± = 0, the ratio A/hil assumes a zero of Jo, i.e., one of the values 
2.405. .,5.520.. ,8.654....... 

So far, we ignored that the driving also affects the coherent dynamics and, thus, the pulse sequence 
of the CNOT operation needs a modification: At the working points of the driven system, the propagator 
becomes U c tf(t,t') — U\(J(t — t')/K) where 

Ufa) = exp(-i<^M ) (42) 

is the propagator corresponding to the ideal Ising interaction 137t . This allows one to implement the 
alternative CNOT operation depicted in Fig.|3p [39,40,59]. Note that the interaction time tj = irh/AJ is 
the same as for the original gate operation in Fig. [3^. Since £/ac(27t/0) is the identity [cf. Eq. J24H . we 
assume for convenience that tj is an integer multiple of the driving period 27r/f2, i.e., ft = 8kJ/H with 
integer k . 

4.3 Numerical solution 

In order to confirm our analytical results, we compute the dissipative propagator as described in Section 
12.31 Thereby, we restrict ourselves to purely harmonic driving f(t) = Acos(ttt)/2. The resulting purity 
loss during the interaction time tj is depicted in Fig. [5] We find that for large temperatures, k-gT > J, 
decoherence is fairly independent of the driving. This behavior changes as the temperature is lowered: 
Once k?,T < J, the purity loss is significantly reduced whenever the ratio A/Kl is close to a zero of the 
Bessel function Jo. Both observations confirm the preceding analytical estimates. The behavior at the first 
working point A ps 2.4 Kl is depicted in Fig. [(J. For relatively low driving frequencies, we find the purity 
loss being proportional to J/O. This significant deviation from the analytical result for small £1 relates to 
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Fig. 5 Purity loss during the interaction time tj as a function of the driving amplitude. The driving 
frequency is 51 = 32 J jh, and the dissipation strength 27ra = 0.01. For A = 0, the undriven situation is 
reproduced. 
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Fig. 6 (a) Purity loss (1 — V) for a pair of qubits with Heisenberg interaction as a function of the driving 
frequency for the temperatures ksT — 0.1J (dashed), 0.01J (solid), and 0.001J (dash-dotted). The dotted 
lines mark the analytical estimate 1 — V(tj) « Tisingt j. (b) Corresponding fidelity defect 1 — T. 



the fact that the low-frequency regime is not within the scope of our rotating-wave approximation which 
assumes ft to be the largest frequency scale. With increasing driving frequency, the discrepancy decreases 
until finally decoherence is dominated by thermal noise oc T and the numerical solution confirms the 
analytical results. Figure 0reveals that the accuracy of our analytical estimates increases with the driving 
frequency: While for the relatively low frequency 51 = 16 J/h, the purity loss is still close to the one of the 
undriven Heisenberg gate, it converges in the limit 51 — > oo to the value obtained for Ising interaction (cf. 
the dashed line in Fig.|4j. 

Still, there remains one caveat: The gate operation in Fig. |3j> relies on the fact that the static effective 
Hamiltonian -ff qu bits describes the dynamics of the driven system sufficiently well — any discrepancy re- 
sults in a coherent error. Therefore, we still have to justify that such coherent errors are sufficiently small. 
As a measure, we employ the so-called fidelity [47] 



T = tr[pidealp(*j)], 



(43) 



S. Kohler and P. Hanggi: Improving the purity of one- and two-qubit gates by AC fields 



13 



10' 




10 



I 



















0.001 



0.01 



0.1 



1 



10 



100 



k B T/J 



Fig. 7 Purity loss shown in Fig.[6]as a function of the temperature. The driving amplitude A ~ 2.4 hQ, is 
adjusted such that 1 — V assumes its first minimum; cf. Fig. [5] All other parameters are as in Fig.|6| 



which is defined as the overlap between the real outcome of the operation, p(tj), and the desired final 
state pidcai = Ui(tt / 4) pinll] ■ (7r/4) in the average over all pure initial states. The time-evolution (1421 of 
the ideal Ising qubit-qubit interaction is characterized by T = 1. Figure |6j) demonstrates that the fidelity 
defect 1 — T at the first working point is even smaller than the purity loss. Thus, we can conclude that 
coherent errors are not of a hindrance. 

4.4 Implementation with quantum dots 

Figure0indicates that the benefits of a "preferably infinitely" large driving frequency and amplitude. Thus, 
a crucial question is whether sufficiently large values are experimentally within reach. For spin qubits in 
quantum dots [24] a typical exchange coupling is J = 0.1 meV which for a temperature T = 10 mK 
corresponds to the solid lines in Figs. [5] and |6] These results demonstrate that driving with the feasible 
frequency £1 = 2tt x 100 J/H « 10 12 Hz and amplitude A = lOmeV already reduces the purity loss by 
two orders of magnitude while the fidelity loss stays at a tolerable level. 

5 Conclusions 

We have investigated the influence of oscillating fields on the coherence properties of one- and two-qubit 
gate operations for three different cases for which it is beneficial. The first case constitutes a continuous- 
wave version of dynamical decoupling of a single qubit from its environment. It is characterized by a 
driving Hamiltonian that does not commute with the bath coupling. There, we have found that a low- 
frequency driving is rather destructive because it generally even enforces decoherence. However, once 
the frequency exceeds the bath cutoff, the coherence properties recover and are finally significantly im- 
proved, especially at high temperatures. Since such a dynamical decoupling by a harmonic driving allows 
higher driving frequencies than the pulsed version, this form of coherence stabilization bears interesting 
perspectives for applications. 

A second possibility for manipulating the decoherence of a single qubit is provided by the physics of 
coherent destruction of tunneling. For such a driving, we have found that the coherence stabilization results 
from the fact that the driving shifts the coherent long-time dynamics of the qubit towards lower frequencies. 
There, the spectral density of an ohmic bath is lower and, consequently, the effective dissipation is weaker. 
This implies that decoherence is generally reduced — most significantly at low temperatures. For high 
temperatures, however, the lower spectral density is counterbalanced by an increasing thermal noise, such 
that in this regime decoherence is essentially not influenced the driving. 
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For two qubits interacting via a Heisenberg exchange coupling, a suited AC field turns the interaction 
into an effective Ising interaction, which is much less sensitive to decoherence. For qubits with such a 
Heisenberg interaction, like e.g. spin qubits, this suggests the following coherence stabilization protocol: 
Use for the CNOT operation a pulse sequence that is suitable for Ising interaction with the latter being 
realized by a Heisenberg interaction with a proper additional AC field. This coherence stabilization scheme 
differs from previous proposals in two respects: First, it is different from dynamical decoupling because 
the driving commutes with the bath coupling. By contrast, the central idea of our scheme is to suppress 
rather the coherent system dynamics "transverse" to this sensitive system coordinate. Consequently, the 
bit-flip noise acts as pure phase noise, which is proportional to the temperature. Cooling, thus, enables 
a further coherence gain. The second difference is that the proposed scheme eliminates also the noise 
stemming from the spectral range above the driving frequency and, thus, is particularly valuable for ohmic 
noise spectra with large cutoff frequencies. Moreover, the driven system still allows one to perform the 
desired CNOT operation with high fidelity and within the same operation time as in the absence of the 
control field. Hence, the gained coherence time fully contributes to the number of feasible gate operations. 
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A Average over all pure states 

In this appendix, we derive formulas for the evaluation of expressions of the type tr(pA) and tr(pApB) 
in an ensemble average over all pure states p = \ijj)(if;\. The state is an element of an A^-dimensional 
Hilbert space. Decomposed into an arbitrary orthonormal basis set {|«)} n =i...jv, it reads 

N 

|^)=X)c n |n>. ( 44 ) 

n=l 

where the only restriction imposed on the coefficients c n is the normalization (ip\ip) = ^ n \ c n\ 2 = 1- 
Hence the ensemble of all pure states is described by the distribution 

9 N N 

P(c u ...,c N ) = -^^ s(l-J2\cn\ 2 )- (45) 

The prefactor on the right-hand side of Eq. d45i has been determined from the normalization 

J d 2 c 1 ...d 2 c N P(c 1 ,...,c N ) = 1 (46) 

of the distribution and j d 2 c denotes integration over the real and the imaginary part of c. We emphasize 
that P(ci, . . . ,cn) is invariant under unitary transformations of the state \tp). The computation of the 
ensemble averages of the coefficients with the distribution J45t is straightforward and yields 

CmC* = —5 mn , (47) 

CmC n C m 'C n / — — _j_ J\ (^ mn ^ m ' n ' ^fnn' $nrn' ) • (48) 
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Using these expressions, we consequently find for the ensemble averages of the expressions tr(pA) and 
ti(pApB) the results 

Ir^-IVplW-^, (49) 

tr(A) tr(B) + tr(AB) 



tv(pApB) = (i,\A\i,){i;\B\i>) = v N{N + 1) "> < 5 °) 

which have been used to obtain the purity decay rate dlOi . 
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